Lattice platforms for performing quantum computations

ABSTRACT

Apparatus and methods for performing quantum computations are disclosed. Such apparatus and methods may include identifying a first quantum state of a lattice having a system of quasi-particles disposed thereon, moving the quasi-particles within the lattice according to at least one predefined rule, identifying a second quantum state of the lattice after the quasi-particles have been moved, and determining a computational result based on the second quantum state of the lattice. Various platforms can be used to physically implement such a quantum computer. Platforms include an optical lattice, a Josephson junction array, a quantum dot, and a crystal structure. Each platform comprises an appropriate array of associated sites that can be used to approximate a desired Kagome geometry. A charge controller is desirably electrically coupled to the platform so that the array may be manipulated as desired. A computing device may be coupled to the charge controller and/or the platform to provide instructions and support, record measurements, and make calculations, for example.

CROSS-REFERENCE TO RELATED APPLICATIONS

The subject matter disclosed and claimed herein is related to the subject matter disclosed and claimed in U.S. patent applications Ser. Nos. 10/909,005, filed Jul. 30, 2004, 10/930,640, filed on even date herewith, and 10/931,082, filed on even date herewith. The disclosure of each of the above-referenced U.S. patent applications is incorporated herein in its entirety.

FIELD OF THE INVENTION

This invention relates in general to the field of quantum computing. More particularly, this invention relates to topological quantum computing.

BACKGROUND OF THE INVENTION

Since the discovery of the fractional quantum Hall effect in 1982, topological phases of electrons have been a subject of great interest. Many abelian topological phases have been discovered in the context of the quantum Hall regime. More recently, high-temperature superconductivity and other complex materials have provided the impetus for further theoretical studies of and experimental searches for abelian topological phases. The types of microscopic models admitting such phases are now better understood. Much less is known about non-abelian topological phases. They are reputed to be obscure and complicated, and there has been little experimental motivation to consider non-abelian topological phases. However, non-abelian topological states would be an attractive milieu for quantum computation.

It has become increasingly clear that if a new generation of computers could be built to exploit quantum mechanical superpositions, enormous technological implications would follow. In particular, solid state physics, chemistry, and medicine would have a powerful new tool, and cryptography also would be revolutionized.

The standard approach to quantum computation is predicated on the quantum bit (“qubit”) model in which one anticipates computing on a local degree of freedom such as a nuclear spin. In a qubit computer, each bit of information is typically encoded in the state of a single particle, such as an electron or photon. This makes the information vulnerable. If a disturbance in the environment changes the state of the particle, the information is lost forever. This is known as decoherence—the loss of the quantum character of the state (i.e., the tendency of the system to become classical). All schemes for controlling decoherence must reach a very demanding and possibly unrealizable accuracy threshold to function.

Topology has been suggested to stabilize quantum information. A topological quantum computer would encode information not in the conventional zeros and ones, but in the configurations of different braids, which are similar to knots but consist of several different threads intertwined around each other. The computer would physically weave braids in space-time, and then nature would take over, carrying out complex calculations very quickly. By encoding information in braids instead of single particles, a topological quantum computer does not require the strenuous isolation of the qubit model and represents a new approach to the problem of decoherence.

In 1997, there were independent proposals by Kitaev and Freedman that quantum computing might be accomplished if the “physical Hilbert space” V of a sufficiently rich TQFT (topological quantum field theory) could be manufactured and manipulated. Hilbert space describes the degrees of freedom in a system. The mathematical construct V would need to be realized as a new and remarkable state for matter and then manipulated at will.

In 2000, Freedman showed that some extraordinarily complicated local Hamiltonian H can be written down whose ground state is V. But this H is an existence theorem only, and is far too complicated to be the starting point for a physical realization.

In 2002, Freedman showed a Hamiltonian involving four-body interactions and stated that after a suitable perturbation, the ground state manifold of H will be the desired state V. This H is less complex than the previously developed H, but it is still only a mathematical construct. One does not see particles, ions, electrons, or any of the prosaic ingredients of the physical world in this prior art model. A Hamiltonian is an energy operator that describes all the possible physical states (eigenstates) of the system and their energy values (eigenvalues).

Freedman further defined the notion of d-isotopy, and showed that if it can be implemented as a ground state of a reasonable Hamiltonian, then this would lead to V and to topological quantum computation. Isotopy is defined as deformation, and two structures that are isotopic are considered to be the same. As shown in the toruses 1 and 2 of FIGS. 1A and 1B, respectively, for example, X and X′ are isotopic, because one may be gradually deformed into the other. In d-isotopy, small circles can be absorbed as a factor=d. Such closed curves are referred to as multicurves or multiloops. Loop X″ in FIG. 1C (winding around torus 3) is not d-isotopic to X or X′. Loops that are unimportant (because, e.g., they comprise a contractible circle) are called trivial loops and it is desirable to remove, as well as count them. Whenever a trivial loop is removed, the picture is multiplied by “d”. In other words, if two multiloops are identical except for the presence of a contractible circle, then their function values differ by a factor of d, a fixed positive real number. It has been shown that d=2 cos π/(k+2), where k is a level such as 1, 2, 3, etc. which is a natural parameter of Cherns-Simons theory.

According to Freedman, the parameter d can take on only the “special” values: 1, root2, golden ratio, root3 . . . 2 cosπ/(k+2) (where k is a natural number). At d=1, the space V becomes something already known, if not observed in solid state physics. For d>1, V is new to the subject. Freedman, et al., later showed that d-isotopy is explicable by field theory and that multiloops as domain walls can be alternately interpreted as Wilson loop operators. Thus, d-isotopy is a mathematical structure that can be imposed on the multiloops, and is based on Cherns-Simons theory.

An exotic form of matter is a fractional quantum Hall fluid. It arises when electrons at the flat interface of two semiconductors are subjected to a powerful magnetic field and cooled to temperatures close to absolute zero. The electrons on the flat surface form a disorganized liquid sea of electrons, and if some extra electrons are added, quasi-particles called anyons emerge. Unlike electrons or protons, anyons can have a charge that is a fraction of a whole number.

The fractional quantum Hall fluids at one-third filling (of the first Landau level) are already a rudimentary (abelian) example of the V of a TQFT. To effect quantum computation, it would be desirable to construct states more stable and more easily manipulated than FQHE (fractional quantum Hall effect) fluids.

One property of anyons is that when they are moved around each other, they remember in a physical sense the knottedness of the paths they followed, regardless of the path's complexity. It is desirable to use anyons in a system with complex enough transformations, called non-abelian transformations, to carry out calculations in a topological quantum computation system.

In view of the foregoing, there is a need for systems and methods that overcome the limitations and drawbacks of the prior art.

SUMMARY OF THE INVENTION

Apparatus and methods according to the invention may include identifying a first quantum state of a lattice having a quasi-particle disposed thereon, moving the quasi-particle within the lattice according to at least one predefined rule, identifying a second quantum state of the lattice after the quasi-particle has been moved, and determining a computational result based on the second quantum state of the lattice. The quasi-particle may be a non-abelian anyon, for example.

Various platforms can be used to physically implement an exemplary quantum computer as set forth above. Exemplary platforms include an optical lattice, a Josephson junction array, a quantum dot, and a crystal structure. Each platform comprises an appropriate array of associated sites that can be used to approximate a desired Kagome geometry. A charge controller is desirably electrically coupled to the platform so that the array may be manipulated as desired. A computing device may be coupled to the charge controller and/or the platform to provide instructions and support, record measurements, and make calculations, for example.

Additional features and advantages of the invention will be made apparent from the following detailed description of illustrative embodiments that proceeds with reference to the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

The foregoing summary, as well as the following detailed description of preferred embodiments, is better understood when read in conjunction with the appended drawings. For the purpose of illustrating the invention, there is shown in the drawings exemplary constructions of the invention; however, the invention is not limited to the specific methods and instrumentalities disclosed. In the drawings:

FIGS. 1A, 1B, and 1C are diagrams useful in describing isotopy;

FIG. 2 is a diagram of an exemplary Kagome lattice in accordance with the present invention;

FIG. 3 is a diagram of an exemplary Kagome lattice in accordance with the present invention;

FIG. 4 is a diagram of an exemplary lattice that is useful for describing aspects of the present invention;

FIG. 5 is a diagram of an exemplary lattice useful for describing dimer moves in accordance with the present invention;

FIG. 6 is a diagram of an exemplary system including an exemplary platform that can be used to physically implement aspects of the present invention; and

FIG. 7 is a block diagram showing an exemplary computing environment in which aspects of the invention may be implemented.

DETAILED DESCRIPTION OF ILLUSTRATIVE EMBODIMENTS

General Description of Quantum Computing

Apparatus and methods according to the invention include identifying a first quantum state of a lattice having a quasi-particle disposed thereon, moving the quasi-particle within the lattice according to at least one predefined rule, identifying a second quantum state of the lattice after the quasi-particle has been moved, and determining a computational result based on the second quantum state of the lattice.

In accordance with the present invention, “realistic” microscopics can be provided for d-isotopy, wherein “realistic” refers to atoms and electrons, for example, as opposed to loops and curves. More particularly, realistic means physical degrees of freedom, with local interaction (the particles are nearby to each other and know of each other's presence). The interactions are potential energy costs of bringing atoms near each other. Existing physical relationships, such as van der waals forces, may be used, along with other physical characteristics, such as tunneling amplitudes. It is desirable to obtain a physical embodiment of the mathematical construct of d-isotopy. In other words, it is desirable to turn the abstract description of the d-isotopy into known physical processes (e.g., tunneling, repulsion, Coulomb interaction).

To go from abstract d-isotopy to real physics, multiloops (i.e., multicurves) are implemented. These multiloops are desirably implemented as broken curves (i.e., “dimer covers”, which is a standard term of art in physics). The following rules are desirably implemented as well: (1) isotopy rule—the amplitude of a picture (or a portion of picture) does not change as it is bent; and (2) d rule—to account for the value associated with a small loop. The rules are turned into “fluctuations”, which are terms in the Hamiltonian that relate to different states (e.g., flickering between two states).

An exemplary embodiment is directed to an extended Hubbard model with a two-dimensional Kagome lattice and a ring-exchange term, described further below. A dimer cover (such as a one-sixth filled Kagome lattice) along with a topology (such as loops) and rules (such as the d-isotopy rule) are desirably comprised within an exemplary extended Hubbard model. Exemplary particles used in conjunction with an exemplary extended model can be bosons or spinless fermions.

At a filling fraction of one-sixth, the model is analyzed in the lowest non-vanishing order of perturbation theory. For a lattice populated with a certain percentage of particles (e.g., electrons), the particles will naturally dissipate to form the “perfect” arrangement. Thus, if one-sixth of the edges of the triangular lattice are filled with particles (e.g., electrons), a perfect matching of the electrons results (e.g., as the result of the Coulomb repulsion).

An exactly soluble point is determined whose ground state manifold is the extensively degenerate “d-isotopy space”, a precondition for a certain type of non-abelian topological order. Near the values d=2 cos π/(k+2), this space collapses to a stable topological phase with anyonic excitations closely related to SU(2) Chern-Simons theory at level k.

A class of Hamiltonians produces a ground state manifold V=d-isotopy (in some cases “weak d-isotopy”) from standard physical processes. The Hamiltonian is an “extended Hubbard model” in certain cases with an additional “ring exchange term”. A parameter domain is described along with a lattice on which an exemplary model operates. An exemplary lattice is the Kagome lattice. An exemplary model as described herein can serve as a blueprint for the construction of phases of matter Vd, which is useful for quantum information processing. The present invention opens the topological path to a quantum computer.

The extended Hubbard model formalizes the Hamiltonian (kinetic+potential energy) as a tunneling term plus energy costs which are dependent on individual particle location (“on sight potential”) and pair-wise locations (e.g., Coulomb repulsion). For this exemplary model, a ring exchange term may be added, which models collective rotations of large groups of particles (e.g., groups of four particles). According to exemplary embodiments, the cases in which the particles are fermions or bosons are treated separately. The particles can be electrons, Cooper pairs, neutral atoms, phonons, or more exotic “electron fractions” such as chargeons, for example. It is noted that spin is not taken into account in the exemplary embodiments, and so if the particle is specialized to be an electron, either spin must be frozen with a magnetic field or a hyperfine splitting of levels tolerated. It is further noted that a formal transformation may be performed such that these occupation based models may be regarded as purely spin models.

The Kagome lattice is well known in condensed matter physics to support highly “frustrated” spin models with rather mysterious ground states. A Kagome lattice is formed using the centers of the edges of the triangular lattice, and resembles a plurality of hexagons with triangles connecting them. An exemplary Kagome lattice used herein is shown with distinguished sublattices, as shown in FIG. 2.

The Hamiltonian given by:

$\begin{matrix} \begin{matrix} {H = {{\sum\limits_{i}{\mu_{i}n_{i}}} + {U_{0}{\sum\limits_{i}\; n_{i}^{2}}} + {U{\sum\limits_{{({i,j})} \in O}{n_{i}n_{j}}}} + {\sum\limits_{{{({i,j})} \in {\vartriangleright \vartriangleleft}},{\notin O}}\;{V_{ij}n_{i}n_{j}}} -}} \\ {{\sum\limits_{({i,j})}{t_{i,j}\left( {{c_{i}^{\phi}c_{j}} + {c_{j}^{\phi}c_{i}}} \right)}} + {{Ring}.}} \end{matrix} & (1) \end{matrix}$ is the occupation number on site i, and μ_(i) is the corresponding chemical potential. U₀ is the usual onsite Hubbard energy U₀ (superfluous for spinless fermions). U is a (positive) Coulomb penalty for having two particles on the same hexagon while V_(ij) represents a penalty for two particles occupying the opposite corners of “bow-ties” (in other words, being next-nearest neighbors on one of the straight lines). Allowing for the possibility of inhomogeneity, not all V_(ij) are assumed equal. Specifically, define v^(c) _(ab)=V_(ij), where a is the color of site (i), b is the color of (j), and c is the color of the site between them. In the lattice of FIG. 2, each of the following is possible distinct, v^(g) _(bb), v^(b) _(bb), V^(g) _(bg), v^(b) _(rb), and v^(b) _(rg), where rεR, gεG, and bεB=K\(R∪G), t_(ij) is the usual nearest-neighbor tunnelling amplitude which is also assumed to depend only on the color of the environment: t_(ij)≡t_(cab) where c now refers to the color of the third site in a triangle. “Ring” is a ring exchange term—an additional kinetic energy term which is added to the Hamiltonian on an ad hoc basis to allow correlated multi-particle hops which “shift” particles along some closed paths.

The onsite Hubbard energy U₀ is considered to be the biggest energy in the problem, and it is set to infinity, thereby restricting the attention to the low-energy manifold with sites either unoccupied or singly-occupied. The rest of the energies satisfy the following relations: U>>t_(ij), V_(ij), μ_(i).

Equations may be derived to the second order in perturbation theory for the ground state manifold of this extended Hubbard model. Such equations are provided below. The solutions describe parameter regimes within the Hubbard model for the existence of the phase called d-isotopy. Technically, the dotted sublattice of the Kagome is desirably altered by a system of defects, or alternatively a “ring exchange term” is desirably introduced to achieve d-isotopy; the bare model yields the ground state manifold “weak d-isotopy”, but this distinction is only technical. Both “weak” and “ordinary” behave similarly at the next step.

It is known from the theory of C*-algebras that for “special” d=2 cos π/(k+2), d-isotopy has a unique symmetry which if broken will relieve the extensive degeneracy of the ground state manifold to a mere finite degeneracy (which depends on topology and boundary conditions). It is further known that once this symmetry is broken, the result is the topological phase Vd. Vd functions as a universal quantum computer. In accordance with the present invention, under a large domain of perturbations, this symmetry will be broken. Thus, the Hamiltonians in accordance with the present invention are blueprints for d-isotopy and hence, after perturbation, for a universal quantum computational system. Such a system has a degree of topological protection from decoherence, which is an enemy of quantum computation. The strength of this protection will depend on the size of the spectral gap above Vd—a quantity which is difficult to compute but can be bounded from the energy scales of any given implantation. There is no a priori basis (such as exists in the FQHE) for asserting that this spectral gap will be small.

Although the equations first solve for an occupation model with Kagome geometry, the invention extends to any physical implementation of the quantum doubles of SU(2) Chern-Simons theory (and their Kauffman-Turaev-Viro) variants which proceed by breaking the symmetry inherent in d-isotopy at its “special” values. Such implementations may be based on lattice models with spin or occupation degrees of freedom or based on field theory.

The non-abelian topological phases which arise are related to the doubled SU(2) Chern-Simons theories described in the prior art. These phases are characterized by (k+1)²-fold ground state degeneracy on the torus T2 and should be viewed as a natural family containing the topological (deconfined) phase of Z₂ gauge theory as its initial element, k=1. For k≧2, the excitations are non-abelian. For k=3 and k≧5, the excitations are computationally universal.

The conditions that a microscopic model should satisfy in order to be in such a topological phase are described. It is useful to think of such a microscopic model as a lattice regularization of a continuum model whose low energy Hilbert space may be described as a quantum loop gas. More precisely, a state is defined as a collection of non-intersecting loops. A Hamiltonian acting on such state can do the following: (i) the loops can be continuously deformed—this “move” is referred to as an isotopy move; (ii) a small loop can be created or annihilated—the combined effect of this move and the isotopy move is referred to as ‘d-isotopy’; and (iii) when exactly k+1 strands come together in some local neighborhood, the Hamiltonian can cut them and reconnect the resulting “loose ends” pairwise so that the newly-formed loops are still non-intersecting.

More specifically, in order for this model to be in a topological phase, the ground state of this Hamiltonian should be a superposition of all such pictures with the additional requirements that (i) if two pictures can be continuously deformed into each other, they enter the ground state superposition with the same weight; (ii) the amplitude of a picture with an additional loop is d times that of a picture without such loop; and (iii) this superposition is annihilated by the application of the Jones-Wenzl (JW) projector that acts locally by reconnecting k+1 strands. It should be noted that, as described herein, the particular form of these projectors is highly constrained and leads to a non-trivial Hilbert space only for special values of d=±2 cos π/(k+2). A Hamiltonian is constructed which enforces d-isotopy for its ground state manifold (GSM).

An exemplary model is defined on the Kagome lattice shown in FIG. 3, which is similar to that shown in FIG. 2. The sites of the lattice are not completely equivalent, and two sublattices are shown, as represented by R (red) and G (green). In FIG. 3, solid dots and dashed lines represent sites and bonds of the Kagome lattice K with the special sublattices R and G. Solid lines define the surrounding triangular lattice.

A dimer cover refers to every vertex meeting exactly one (i.e., one and only one) edge. Dimer covers arise physically (naturally) by repulsion. For example, electrons repel each other to form dimer covers.

As shown, green is a perfect match (as defined above), and red is a second perfect match. Two perfect matches result in multiloops that alternate in color (green, red, green, red, etc.).

It is possible that a green dimer and a red dimer cover the same edge. A solution is to consider it a very short loop of length two (a red and green—assume one of the red or green is slightly displaced so it is a “flat” loop). The multiloops will have alternated these dimers.

Encoding is desirable for an exemplary quantum computing model. The loops are encoded in dimers, and the rules are encoded in particle (e.g., electronic) interactions, such as repulsion or tunneling, for example. In certain situations, when electrons get too close, the system passes from the ground state to the excited state and back to the ground state. This is a virtual process seen in second order perturbation theory and is used to build fluctuations; it generates a deformation.

The “undoped” system corresponds to the filling fraction one-sixth (i.e., Np≡Σi ni=N/6, where N is the number of sites in the lattice). The lowest-energy band then consists of configurations in which there is exactly one particle per hexagon, hence all U-terms are set to zero. These states are easier to visualize if a triangular lattice T is considered whose sites coincide with the centers of hexagons of K, where K is a surrounding lattice for T. Then a particle on K is represented by a dimer on T connecting the centers of two adjacent hexagons of K.

The condition of one particle per hexagon translates into the requirement that no dimers share a site. In the ⅙-filled case, this low-energy manifold coincides with the set of all dimer coverings (perfect matchings) of T. The “red” bonds of T (the ones corresponding to the sites of sublattice R) themselves form one such dimer covering, a so-called “staggered configuration”. This particular covering is special: it contains no “flippable plaquettes”, or rhombi with two opposing sides occupied by dimers. See FIG. 4, which shows a triangular lattice T obtained from K by connecting the centers of adjacent hexagons. The bonds corresponding to the special sublattices R and G are shown in dashed and dotted lines, respectively. Triangles with one red side are shaded.

Therefore, particles live on bonds of the triangular lattice and are represented as dimers. In particular, a particle hop corresponds to a dimer “pivoting” by 60 degrees around one of its endpoints. V_(ij)=V^(c) _(ab) is now a potential energy of two parallel dimers on two opposite sides of a rhombus, with c being the color of its short diagonal.

Desirably, the triangular lattice is not bipartite. On the edges of a bipartite lattice, the models will have an additional, undesired, conserved quantity (integral winding numbers, which are inconsistent with the JW projectors for k>2), so the edge of the triangular lattice gives a simple realization.

Because a single tunneling event in D leads to dimer “collisions” (two dimers sharing an endpoint) with energy penalty U, the lowest order at which the tunneling processes contribute to the effective low-energy Hamiltonian is 2. At this order, the tunneling term leads to two-dimer “plaquette flips” as well as renormalization of bare onsite potentials due to dimers pivoting out of their positions and back.

By fixing R as in FIG. 3, without small rhombi with two opposite sides red, as the preferred background dimerization, the fewest equations are obtained along with ergodicity under a small set of moves. Unlike in the usual case, the background dimerization R is not merely a guide for the eyes, it is physically distinguished: the chemical potentials and tunneling amplitudes are different for bonds of different color.

Exemplary elementary dimer moves that preserve the proper dimer covering condition include plaquette flips, triangle moves, and bow-tie moves. A plaquette (rhombus) flip is a two-dimer move around a rhombus made of two lattice triangles. Depending on whether a “red” bond forms a side of such a rhombus, its diagonal, or is not found there at all, the plaquettes are referred to, respectively, as type 1 (or 1′), 2, or 3 (see the lattice diagram of FIG. 5). FIG. 5 shows an overlap of a dimer covering of T (shown in thick black line) with the red covering shown in dashed line) corresponding to the special sublattice R. Shaded plaquettes correspond to various dimer moves described herein.

The distinction between plaquettes of type 1 and 1′ is purely directional: diagonal bonds in plaquettes of type 1 are horizontal, and for type 1′ they are not. This distinction is desirable because the Hamiltonian breaks the rotational symmetry of a triangular (or Kagome) lattice. A triangle move is a three-dimer move around a triangle made of four elementary triangles. One such “flippable” triangle is labeled 4 in FIG. 5. A bow-tie move is a four-dimer move around a “bow-tie” made of six elementary triangles. One such “flippable” bow tie is labeled 5 in FIG. 5.

To make each of the above moves possible, the actual dimers and unoccupied bonds desirably alternate around a corresponding shape. For both triangle and bow-tie moves, the cases when the maximal possible number of “red” bonds participate in their making (2 and 4 respectively) are depicted. Note that there are no alternating red/black rings of fewer than 8 lattice bonds (occupied by at most 4 non-colliding dimers). Ring moves only occur when red and black dimers alternate; the triangle labeled 4 in FIG. 5 does not have a ring term associated with it, but the bow-tie labeled 5 does.

The correspondence between the previous smooth discussion and rhombus flips relating dimerizations of J of T is now described. The surface is now a planar domain with, possibly, periodic boundary conditions (e.g., a torus). A collection of loops is generated by R∪J, with the convention that the dimers of R∩J be considered as length 2 loops or bigons). Regarding isotopy, move 2 is an isotopy from R∪J to R∪J′ but by itself, it does almost nothing. It is impossible to build up large moves from type 2 alone. So it is a peculiarity of the rhombus flips that there is no good analog of isotopy alone but instead go directly to d-isotopy. The following relations associated with moves of type 5 and 1 (1′) are imposed: d ³Ψ

−Ψ

=0_(:)  (2a) dΨ

−Ψ

=0  (2b) because from one to four loops in Equation (2a) is passed, and zero to one loop is passed in Equation (2b).

Having stated the goal, the effective Hamiltonian H: D→ D

-   -   is derived on the span of dimerizations. The derivation is         perturbative to the second order in ε where ε=t^(r)         _(bb)/U=t^(b) _(gb)/U. Additionally, t^(b) _(rb)/U=c₀ε where c₀         is a positive constant, while t^(g) _(bb)=o(ε) and can be         neglected in the second-order calculations. (In the absence of a         magnetic field, all t's can be made real and hence symmetric         with respect to their lower indices. Also, U=1 for notational         convenience.) Account for all second-order processes, i.e.,         those processes that take us out of D and then back to D. These         amount to off-diagonal (hopping) processes—“plaquette flips” or         “rhombus moves”—as well as diagonal ones (potential energy) in         which a dimer pivots out and then back into its original         position. The latter processes lead to renormalization of the         bare onsite potentials μ_(i), which are adjusted so that all         renormalized potentials are equal up to corrections 0(ε³). The         non-constant part of the effective Hamiltonian comes from the         former processes and can be written in the form:

$\overset{\sim}{H} = {\sum\limits_{I,J}{\left( {{\overset{\sim}{H}}_{I,J} \otimes {II}} \right)\Delta_{I,J}}}$ where ${\overset{\sim}{H}}_{I,J}$

-   -    is a 2×2 matrix corresponding to a dimer move in the         two-dimensional basis of dimer configurations connected by this         move. Δ_(IJ)=1 if the dimerizations I, JεD are connected by an         allowed move, and Δ_(IJ)=0 otherwise.

Therefore, it suffices to specify these 2×2 matrices for the off-diagonal processes. For moves of types (1)-(3), they are given below:

$\begin{matrix} {{{\overset{\sim}{H}}^{(1)} = {\begin{pmatrix} v_{gb}^{b} & {{- 2}t_{rb}^{b}t_{gb}^{b}} \\ {{- 2}t_{rb}^{b}t_{gb}^{b}} & v_{rb}^{b} \end{pmatrix} = \begin{pmatrix} v_{gb}^{b} & {{- 2}c_{0}\varepsilon^{2}} \\ {{- 2}c_{0}\varepsilon^{2}} & v_{rb}^{b} \end{pmatrix}}},} & \left( {3a} \right) \\ {{{\overset{\sim}{H}}^{(1^{\prime})} = {\begin{pmatrix} v_{bb}^{b} & {{- 2}t_{rb}^{b}t_{gb}^{b}} \\ {{- 2}t_{rb}^{b}t_{gb}^{b}} & v_{rb}^{b} \end{pmatrix} = \begin{pmatrix} v_{bb}^{b} & {{- 2}c_{0}\varepsilon^{2}} \\ {{- 2}c_{0}\varepsilon^{2}} & v_{rb}^{b} \end{pmatrix}}},} & \left( {3b} \right) \\ {{{\overset{\sim}{H}}^{(2)} = {\begin{pmatrix} v_{bb}^{b} & {{- 2}\left( t_{bb}^{r} \right)^{2}} \\ {{- 2}\left( t_{bb}^{r} \right)^{2}} & v_{bb}^{r} \end{pmatrix} = \begin{pmatrix} v_{bb}^{r} & {{- 2}\varepsilon^{2}} \\ {{- 2}\varepsilon^{2}} & v_{bb}^{r} \end{pmatrix}}},} & \left( {3c} \right) \\ {{\overset{\sim}{H}}^{(3)} = {\begin{pmatrix} v_{bb}^{g} & {{- 2}\left( t_{bb}^{g} \right)^{2}} \\ {{- 2}\left( t_{bb}^{g} \right)^{2}} & v_{bb}^{g} \end{pmatrix} = {\begin{pmatrix} v_{bb}^{g} & 0 \\ 0 & v_{rb}^{g} \end{pmatrix}.}}} & \left( {3d} \right) \end{matrix}$

H can now be tuned to the “small loop” value d.

${\overset{\sim}{H}}^{(1)} = {{\overset{\sim}{H}}^{(1^{\prime})} \propto \begin{pmatrix} d & {- 1} \\ {- 1} & d^{- 1} \end{pmatrix}}$ is required because these moves change the number of small loops by one. Because a move of type 2 is an isotopy move,

${\overset{\sim}{H}}^{(2)} \propto \begin{pmatrix} 1 & {- 1} \\ {- 1} & 1 \end{pmatrix}$

H ⁽³⁾=0 provided d>1, because it represents a “surgery” on two strands not allowed for k>1. For k=1, on the other hand,

${\overset{\sim}{H}}^{(3)} \propto {\begin{pmatrix} 1 & {- 1} \\ {- 1} & 1 \end{pmatrix}.}$

At level k=1 configurations which differ by such a surgery should have equal coefficients in any ground state vector Ψ while at levels k>1 no such relation should be imposed. Thus, for k>1, the matrix relations (3 a-3 d) yield equations in the model parameters: Types (1)&(1′): v _(gb) ^(b) =v _(bb) ^(b)=2dc ₀ε²  (4a) and v _(rb) ^(b) =v _(rg) ^(b)=2d ⁻¹ c ₀ε²  (4b) Types (2)&(3): v _(bb) ^(r)=2ε² and v _(bb) ^(g)=0  (4c)

Suppose that the Hamiltonian has a bare ring exchange term, ring in Eq. (1):

${\overset{\sim}{H}}^{({Ring})} \propto {\begin{pmatrix} x & {{- c_{3}}\varepsilon^{2}} \\ {{- c_{3}}\varepsilon^{2}} & y \end{pmatrix}.}$ for some constants c3, x, y>0, and consider the additional equations which come from considering Ring as a fluctuation between one loop of length 8 (type 5 in FIG. 5) and four bigons. It follows from equation (2a) that

$\begin{matrix} {{{\overset{\sim}{H}}^{({Ring})} \propto \begin{pmatrix} d^{3} & {- 1} \\ {- 1} & d^{- 3} \end{pmatrix}},{{so}\text{:}}} & \; \\ {{x = {d^{3}c_{3}\varepsilon^{2}}},{y = {d^{- 3}c_{3}{\varepsilon^{2}.}}}} & (5) \end{matrix}$ are the additional equations (beyond equations (4)) to place the model at the soluble point characterized by d-isotopy. It is clear from equation (5) that the diagonal entries of

H ^((Ring))x≠y—perhaps not the most natural choice (although these entries may be influenced by the local chemical environment and thus do not need to be equal). Another possibility exploits the ambiguity of whether a bigon should be considered a loop or not—which allows one to choose x=y.

This construction shows how an extended Hubbard model with the additional ring exchange term (or the equivalent quantum dimer model) can be tuned to have a d-isotopy space as its GSM. Further tuning to its special values then allows the JW projectors to reduce this manifold to the correct ground state corresponding to the unique topological phase associated with this value of d. A simple candidate for a universal quantum computer would be tuned to d=(1+sqrt 5)/2.

Although the exemplary embodiments are described with respect to a triangular lattice, it is contemplated that other lattices, such as a square or hexagonal lattice, may be used in accordance with the present invention. Moreover, the use of irregular lattices is contemplated.

Example Lattice Platforms For Performing Quantum Computations

Various platforms can be used to physically implement a quantum computer as set forth above. FIG. 6 is a diagram of an exemplary system including an exemplary platform 60 that can be used to physically implement aspects of the present invention. Exemplary platforms 60 include an optical lattice, a Josephson junction array, a quantum dot array, and a crystal structure. Each platform 60 comprises an appropriate array 63 of associated sites that can be used to approximate a desired Kagome geometry. A charge controller 65 is desirably electrically coupled to the platform 60 so that the array 63 may be manipulated as desired. A computing device 110, such as that described with respect to FIG. 7, may be coupled to the charge controller and/or the platform to provide instructions and support, record measurements, and make calculations, for example.

Because arriving near an exactly soluble point requires fairly precise control of model parameters, synthetic systems are preferred for a realization of the topological phase Vd. An optical lattice is an example of a synthetic system in which a quantum computer can be implemented. An optical lattice can hold thousands of neutral atoms arranged in approximate Kagome geometry with tunneling amplitudes, and one and two body interactions determined by the optical intensity profile. The optical intensity profile can be tuned to achieve Vd.

In such a system, for example, a laser may be shined through diffraction gratings to get an intensity pattern, such as a line or plane or a three-dimensional pattern. The gratings are chosen to produce a certain pattern. This pattern is provided, for example, as the array 63 on the platform 60.

According to an exemplary embodiment, a two-dimensional pattern (on a plane or platform 60) is obtained by shining a laser through diffraction gratings. The Kagome lattice is desirably disposed on the plane, such that the dark spots are at the points of the lattice. Neutral atoms (e.g., sodium, potassium, etc.) are desirably placed into the plane at a point where they cannot fall any further, using well known techniques. For example, a magnetic trap can be used to position each atom, and then the magnetic field is turned off or otherwise removed. Depending on the type of atom used, it can have an attractive or repulsive force.

Desirably, one-sixth as many atoms as possible sites are placed in the plane (i.e., the filling fraction of the array 63 is one-sixth). The atoms then spread out because of their repulsive forces, and form some dimer covering. The loops of the quantum computer, described above, have been translated into dimers in this exemplary embodiment.

The rules described above are translated into the optical lattice. The rules are related to the repulsion and to tunneling. With respect to tunneling, the potential can be lowered between sites to encourage tunneling (e.g., lower the intensity of light lowers the potential). To lower the repulsion between the pairs, the geometry of the lattice can be distorted (by controlling the light intensity) to make them closer or farther apart, or deepen or make more shallow a pit where the atoms rest. A laser may be used as the charge controller 65 in such an exemplary embodiment.

A Josephson junction array is another type of synthetic system in which a quantum computer can be implemented, to support the phase Vd. A Josephson junction is a type of electronic circuit capable of switching at very high speeds when operated at temperatures approaching absolute zero. A Josephson junction exploits the phenomenon of superconductivity, which is the ability of certain materials to conduct electric current with practically zero resistance. A Josephson junction comprises two superconductors, separated by a nonsuperconducting layer so thin that electrons can cross through the insulating barrier. The flow of current between the superconductors in the absence of an applied voltage is called a Josephson current, and the movement of electrons across the barrier is known as Josephson tunneling.

An array 63 of Josephson junctions can be arranged on a platform 60 and can be made into a Kagome lattice by etching, for example. The rules set forth above are implemented by tuning the structure (e.g., adjusting its shape) using the controller 65, for example. In a Josephson junction array, the Kagome lattice is made of superconducting islands that correspond to the sites of the lattice. The shape of the islands can be controlled by controlling the islands' growth (e.g., using molecular beam epitaxy). In a Josephson junction array, the particles are desirably Cooper pairs.

Another type of synthetic system that can be used in a topological quantum computer is a quantum dot array. A quantum dot is a particle of matter so small that the addition or removal of an electron changes its properties. More particularly, a quantum dot is a confined droplet of electrons, usually on a scale of a few nanometers in diameter. The position of a single electron in a quantum dot might attain several states, so that a quantum dot could represent a byte of data.

An array 63 of quantum dots can be arranged on a platform 60 and can be made into a Kagome lattice by etching, for example. The rules can be implemented by controlling the charge on the array, using a charge controller 65, for example. An exemplary technique for controlling the charge is to implement a capacitor layer above the lattice. A single electron is used. The spin of a single electron can influence the system (it can be considered a perturbation).

Another desirable realization of Vd is as a two-dimensional crystal (e.g., grown on a substrate), having an approximate Kagome geometry. A crystal structure is an example of a natural system in which a quantum computer can be implemented. A natural system is generally more rigid than a synthetic system. A two-dimensional crystal can be grown that has a Kagome geometry. In such a case, the crystal itself would be the platform 60 having a structure in the form of an array 63. One atom can be substituted for another atom. The crystal structure can be squeezed by, for example, cooling or warming the structure, or by subjecting the structure to a magnetic field, using a charge controller 65, for example. The spin of a single electron can influence the system (it can be considered a perturbation). A crystal is not easily tunable, but it may not have to be tuned, if the tolerable region of error (i.e., the tolerable amount of variability) is large enough. An exemplary crystal with a desirable configuration is chromium sarocite.

Exemplary Computing Environment

FIG. 7 illustrates an example of a suitable computing system environment 100 in which the invention may be implemented. The computing system environment 100 is only one example of a suitable computing environment and is not intended to suggest any limitation as to the scope of use or functionality of the invention. Neither should the computing environment 100 be interpreted as having any dependency or requirement relating to any one or combination of components illustrated in the exemplary operating environment 100.

The invention is operational with numerous other general purpose or special purpose computing system environments or configurations. Examples of well known computing systems, environments, and/or configurations that may be suitable for use with the invention include, but are not limited to, personal computers, server computers, hand-held or laptop devices, multiprocessor systems, microprocessor-based systems, set top boxes, programmable consumer electronics, network PCs, minicomputers, mainframe computers, distributed computing environments that include any of the above systems or devices, and the like.

The invention may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. The invention may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network or other data transmission medium. In a distributed computing environment, program modules and other data may be located in both local and remote computer storage media including memory storage devices.

With reference to FIG. 7, an exemplary system for implementing the invention includes a general purpose computing device in the form of a computer 110. Components of computer 110 may include, but are not limited to, a processing unit 120, a system memory 130, and a system bus 121 that couples various system components including the system memory to the processing unit 120. The system bus 121 may be any of several types of bus structures including a memory bus or memory controller, a peripheral bus, and a local bus using any of a variety of bus architectures. By way of example, and not limitation, such architectures include Industry Standard Architecture (ISA) bus, Micro Channel Architecture (MCA) bus, Enhanced ISA (EISA) bus, Video Electronics Standards Association (VESA) local bus, and Peripheral Component Interconnect (PCI) bus (also known as Mezzanine bus).

Computer 110 typically includes a variety of computer readable media. Computer readable media can be any available media that can be accessed by computer 110 and includes both volatile and non-volatile media, removable and non-removable media. By way of example, and not limitation, computer readable media may comprise computer storage media and communication media. Computer storage media includes both volatile and non-volatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data. Computer storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, digital versatile disks (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can accessed by computer 110. Communication media typically embodies computer readable instructions, data structures, program modules or other data in a modulated data signal such as a carrier wave or other transport mechanism and includes any information delivery media. The term “modulated data signal” means a signal that has one or more of its characteristics set or changed in such a manner as to encode information in the signal. By way of example, and not limitation, communication media includes wired media such as a wired network or direct-wired connection, and wireless media such as acoustic, RF, infrared and other wireless media. Combinations of any of the above should also be included within the scope of computer readable media.

The system memory 130 includes computer storage media in the form of volatile and/or non-volatile memory such as ROM 131 and RAM 132. A basic input/output system 133 (BIOS), containing the basic routines that help to transfer information between elements within computer 110, such as during start-up, is typically stored in ROM 131. RAM 132 typically contains data and/or program modules that are immediately accessible to and/or presently being operated on by processing unit 120. By way of example, and not limitation, FIG. 7 illustrates operating system 134, application programs 135, other program modules 136, and program data 137.

The computer 110 may also include other removable/non-removable, volatile/non-volatile computer storage media. By way of example only, FIG. 7 illustrates a hard disk drive 140 that reads from or writes to non-removable, non-volatile magnetic media, a magnetic disk drive 151 that reads from or writes to a removable, non-volatile magnetic disk 152, and an optical disk drive 155 that reads from or writes to a removable, non-volatile optical disk 156, such as a CD-ROM or other optical media. Other removable/non-removable, volatile/non-volatile computer storage media that can be used in the exemplary operating environment include, but are not limited to, magnetic tape cassettes, flash memory cards, digital versatile disks, digital video tape, solid state RAM, solid state ROM, and the like. The hard disk drive 141 is typically connected to the system bus 121 through a non-removable memory interface such as interface 140, and magnetic disk drive 151 and optical disk drive 155 are typically connected to the system bus 121 by a removable memory interface, such as interface 150.

The drives and their associated computer storage media provide storage of computer readable instructions, data structures, program modules and other data for the computer 110. In FIG. 7, for example, hard disk drive 141 is illustrated as storing operating system 144, application programs 145, other program modules 146, and program data 147. Note that these components can either be the same as or different from operating system 134, application programs 135, other program modules 136, and program data 137. Operating system 144, application programs 145, other program modules 146, and program data 147 are given different numbers here to illustrate that, at a minimum, they are different copies. A user may enter commands and information into the computer 110 through input devices such as a keyboard 162 and pointing device 161, commonly referred to as a mouse, trackball or touch pad. Other input devices (not shown) may include a microphone, joystick, game pad, satellite dish, scanner, or the like. These and other input devices are often connected to the processing unit 120 through a user input interface 160 that is coupled to the system bus, but may be connected by other interface and bus structures, such as a parallel port, game port or a universal serial bus (USB). A monitor 191 or other type of display device is also connected to the system bus 121 via an interface, such as a video interface 190. In addition to the monitor, computers may also include other peripheral output devices such as speakers 197 and printer 196, which may be connected through an output peripheral interface 195.

The computer 110 may operate in a networked environment using logical connections to one or more remote computers, such as a remote computer 180. The remote computer 180 may be a personal computer, a server, a router, a network PC, a peer device or other common network node, and typically includes many or all of the elements described above relative to the computer 110, although only a memory storage device 181 has been illustrated in FIG. 7. The logical connections depicted include a LAN 171 and a WAN 173, but may also include other networks. Such networking environments are commonplace in offices, enterprise-wide computer networks, intranets and the internet.

When used in a LAN networking environment, the computer 110 is connected to the LAN 171 through a network interface or adapter 170. When used in a WAN networking environment, the computer 110 typically includes a modem 172 or other means for establishing communications over the WAN 173, such as the internet. The modem 172, which may be internal or external, may be connected to the system bus 121 via the user input interface 160, or other appropriate mechanism. In a networked environment, program modules depicted relative to the computer 110, or portions thereof, may be stored in the remote memory storage device. By way of example, and not limitation, FIG. 7 illustrates remote application programs 185 as residing on memory device 181. It will be appreciated that the network connections shown are exemplary and other means of establishing a communications link between the computers may be used.

As mentioned above, while exemplary embodiments of the present invention have been described in connection with various computing devices, the underlying concepts may be applied to any computing device or system.

The various techniques described herein may be implemented in connection with hardware or software or, where appropriate, with a combination of both. Thus, the methods and apparatus of the present invention, or certain aspects or portions thereof, may take the form of program code (i.e., instructions) embodied in tangible media, such as floppy diskettes, CD-ROMs, hard drives, or any other machine-readable storage medium, wherein, when the program code is loaded into and executed by a machine, such as a computer, the machine becomes an apparatus for practicing the invention. In the case of program code execution on programmable computers, the computing device will generally include a processor, a storage medium readable by the processor (including volatile and non-volatile memory and/or storage elements), at least one input device, and at least one output device. The program(s) can be implemented in assembly or machine language, if desired. In any case, the language may be a compiled or interpreted language, and combined with hardware implementations.

The methods and apparatus of the present invention may also be practiced via communications embodied in the form of program code that is transmitted over some transmission medium, such as over electrical wiring or cabling, through fiber optics, or via any other form of transmission, wherein, when the program code is received and loaded into and executed by a machine, such as an EPROM, a gate array, a programmable logic device (PLD), a client computer, or the like, the machine becomes an apparatus for practicing the invention. When implemented on a general-purpose processor, the program code combines with the processor to provide a unique apparatus that operates to invoke the functionality of the present invention. Additionally, any storage techniques used in connection with the present invention may invariably be a combination of hardware and software.

While the present invention has been described in connection with the preferred embodiments of the various figures, it is to be understood that other similar embodiments may be used or modifications and additions may be made to the described embodiments for performing the same functions of the present invention without deviating therefrom. Therefore, the present invention should not be limited to any single embodiment, but rather should be construed in breadth and scope in accordance with the appended claims. 

1. A method for performing quantum computation, comprising: disposing a lattice array of non-Abelian anyons on a platform, the lattice array comprising a system of real particles, the non-Abelian anyons being excitations of least energy states of the system of real particles that form the lattice array; identifying a first quantum state of the lattice array; moving the non-Abelian anyons within a 2D+1 space-time of the lattice array according to at least one predefined rule; identifying a second quantum state of the lattice array after the non-Abelian anyons have been moved; and determining a computational result based on the second quantum state of the lattice array.
 2. The method of claim 1, wherein the platform is an optical lattice.
 3. The method of claim 2, wherein moving the non-Abelian anyons within the lattice array comprises affecting a geometric relationship among the real particles that form the lattice array.
 4. The method of claim 2, wherein disposing the lattice array on the platform comprises shining a laser light through a diffraction grating onto the platform.
 5. The method of claim 2, wherein moving the non-Abelian anyons comprises adjusting the intensity of a light that is shining on the platform.
 6. The method of claim 1, wherein the platform is a Josephson junction array.
 7. The method of claim 6, wherein disposing the lattice array on the platform comprises etching the Josephson junction array.
 8. The method of claim 6, wherein the lattice array comprises superconducting islands disposed on the platform.
 9. The method of claim 8, wherein moving the non-Abelian anyons comprises adjusting the flow of current between the superconducting islands.
 10. The method of claim 1, wherein the platform is a quantum dot array.
 11. The method of claim 10, wherein disposing the lattice array on the platform comprises etching the quantum dot array.
 12. The method of claim 10, wherein moving the non-Abelian anyons comprises controlling a charge on the lattice array.
 13. The method of claim 1, wherein the platform is a crystal structure.
 14. The method of claim 13, wherein disposing the lattice array on the platform comprises growing a two-dimensional crystal having a Kagome geometry.
 15. The method of claim 13, wherein moving the non-Abelian anyons comprises adjusting the temperature of the lattice array.
 16. The method of claim 13, wherein moving the non-Abelian anyons comprises subjecting the lattice array to a magnetic field.
 17. The method of claim 1, wherein the lattice array is a Kagome lattice.
 18. A quantum computing system, comprising: a platform comprising a lattice array having a plurality of non-Abelian anyons disposed thereon, the lattice array comprising a system of real particles, the non-Abelian anyons being excitations of least energy states of the system of real particles that form the lattice array; means for identifying a first quantum state of the lattice array; means for moving the non-Abelian anyons within a 2D+1 space-time of the lattice array according to a set of predefined rules; means for identifying a second quantum state of the lattice array after the non-Abelian anyons are moved; and means for determining a computational result based on the second quantum state of the lattice array.
 19. The system of claim 18, wherein the platform is an optical lattice.
 20. The system of claim 19, wherein the means for moving the non-Abelian anyons within the lattice array comprises means for affecting a geometric relationship among the real particles that form the lattice array.
 21. The system of claim 19, further comprising a laser light and diffraction grating for disposing the lattice array on the platform.
 22. The system of claim 19, wherein the means for moving the non-Abelian anyons comprises a laser light that is shining on the platform.
 23. The system of claim 18, wherein the platform is a Josephson junction array.
 24. The system of claim 23, wherein the lattice array comprises superconducting islands disposed on the platform.
 25. The system of claim 24, wherein the means for moving the non-Abelian anyons comprises means for adjusting the flow of current between the superconducting islands.
 26. The system of claim 18, wherein the platform is a quantum dot array.
 27. The system of claim 26, wherein the means for moving the non-Abelian anyons comprises means for controlling the charge on the lattice array.
 28. The system of claim 18, wherein the platform is a crystal structure.
 29. The system of claim 28, further comprising means for growing a two-dimensional crystal having a Kagome geometry on the platform.
 30. The system of claim 28, wherein the means for moving the non-Abelian anyons comprises one of means for adjusting the temperature of the lattice array and means for subjecting the lattice array to a magnetic field.
 31. The system of claim 18, wherein the lattice array is a Kagome lattice. 